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O |- Abstract. In this paper we consider the problem of simultaneously determin- 

ing the time-dependent thermal diffusivity and the temperature distribution 
in one-dimensional heat equation in the case of nonlocal boundary and inte- 
C _ ) gral overdetermination conditions. We establish conditions for the existence 

f^") , and uniqueness of a classical solution of the problem under considerations. We 

present some results on the numerical solution with an example. 

< 

^ . Introduction 



c3 



Suppose that one need to determine the temperature distribution u(x,t) as 
well as thermal coefficient a(t) simultaneously satisfy the equation 



>■ (0.1) u t = a{t)u xx + F (x,t), < x < 1, < t < T, 

' with the initial condition 

in 
in 

■ the boundary conditions 

o ' 



(0.2) u(x,0) = ip(x), < x< 1, 



(0.3) u(0, t) = u(l, t), u x {\, t)=0, < t < T, 

and the overdetermination condition 

l 

(0.4) J u(x,t)dx = E(t), < t < T. 

o 

The problem of finding the pair {a(i), u(x, t)} in (0.1)-(0.4) will be called an inverse 
problem. 

Denote the domain Qt by 



Q T = {(x, t): < x < 1, < t < T} . 
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Definition 1. The pair{a(t), u(x,t)} from the class C [G^TjxC 2 ' 1 (Qt)^C 1 - (Q t ) 
for which conditions (0.1)- (0.4) are satisfied and a(t) > on the interval [0,T] , is 
called the classical solution of the inverse problem (0.1)-(0.4). 

The parameter identification in a parabolic differential equation from the data 
of integral overdetermination condition plays an important role in engineering and 
physics. ([1,2,3,4,5]) 

Various statements of inverse problems on determination of thermal coefficient 
in one-dimensional heat equation were studied in [4,5,6]. It is important to note 
that in the papers [4,5] the time dependent thermal coefficient is determined by 
nonlocal overdetermination condition's data. Besides, in [1,4] the coefficients of the 
heat equations are determined in the case of nonlocal boundary conditions. 

In the present work, the existence and uniqueness of the classical solution of 
the problem (0.1)-(0.4) is reduced to fixed point principles by applying Fourier 
method. The boundary conditions (0.3) admit the expantions by the system of 
eigenfunctions and associated functions corresponding to the spectral problem. 

The paper organized as follows: 

In Chapter 1, the auxiliary spectral problem which can be obtained by applying 
Fourier method to the problem (0.1)-(0.3) is studied. In Chapter 2, the existence 
of the solution of inverse problem in Qt, and the uniqueness of the solution of the 
inverse problem in Qt (0 < T < T) are shown. Then in Chapter 3, the continuous 
dependence upon the solution of the inverse problem is shown. Finally, in Chapter 
4, the numerical solution for the inverse problem is presented with an example. 

1. The auxiliary spectral problem 
Consider the spectral problem 

X"{x)+\X{x) =0, 0<x< 1, 

(1.1) X(0) = X(1), X'(1)=0. 

This problem is wellknown in [7], as auxiliary spectral problem in solving a 
nonlocal boundary value problem for heat equation by Fourier method. 
It is clear to show that, the problem (1.1) has eigenvalues 

A fe = (27rfc) 2 ,fc = 0, 1,2, ... 

and eigenfunctions 

(1.2) X (x) = 2, X k (x) = 4cos27rfcx, fc = l,2,... 

and the system of functions Xk(x), k = 0, 1, 2, ... is not basis in L2[0, 1]. Complete 
the system Xk(x), k = 0, 1, 2, ... with the associated functions 

(1.3) X k (x) = 4(1 - a;) sin 27rfcx, fc = 1,2,... 

of the problem (1.1). Denote the systems of functions (1.2) and (1.3) as follows: 

(1.4) 

Xq{x) = 2, X 2 k-i(x) — 4cos27rfca;, X 2 k(x) = 4(1 — x) sin27rfcx, fc = l,2, ... . 
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The system of functions Xk(x), k = 0, 1, 2, ... is basis in L2[0, 1]. ([8]) 
The adjoint problem of (1.1) has the form 

Y"(x) + XY(x) = 0, < x < 1, 

(1.5) Y(0) = 0, r'(o)=y'(i). 

Analogously to the system (1.4), the system of eigenfunctions and associated func- 
tions of the problem (1.5) is denoted by 

(1.6) Yo(x) = x, Y2k-i{x) — xcos27rfcx, Y2k(x) = sin27rfcx, fc = l,2, .... 

It is easy to calculate that the systems (1.4) and (1.6) form a biorthonormal 
system on interval [0, 1] , i.e. 

(Xu Y 3 ) = J 1 X l (x)Y 1 (x)dx = % = { J; ■ f • 

2. Existence and Uniqueness of the solution of the inverse problem 

We have the following assumptions on the data of the problem (0.1)-(0.4). 
(Ax) E(t) G C 1 [0,T] , E'(t) < 0, Vt G [0,T] ; 
(A 2 ) tp(x) G C 4 [0, 1] ; 

(1) p(0) = ^(1) = 0, p"(0) = <p"(l),J<p(x)dx = E(0); 

(2) <p 2k > 0,_k = 1,2, 

(A 3 ) F(ar, t)eC (Q T ) ; F(z, t) G C7 4 [0, 1] for arbitrary fixed t £ [0, T] ; 

(1) F(0,t) = F(l,t), F x (l,t) = 0, F xx (0,t) = F xx (l,t); 

T oo T 

(2) F 2 fe(t)>0,fc = 0,1, 2, ....,/ £'(t)dt + £ 4^2fe-2/ F (t)dt>0, 

fc=l 

1 1 
where ip k = J ip(x)Y k (x)dx, F k (t) = J F(x,t)Y k (x)dx, k — 0,1,2,... . 
o o 

Remark 1. There are functions <p, E and F satisfying the assumptions (Ai) — (A3) . 
For example 

(p(x) = (1 — x) sin27ra;, 

E(t) = ^exp(-i), 

F(x,t) = (1 — x) sin27rxexp(3t). 

The main result is presented as follows. 

Theorem 1. Let the assumptions (A\) — (A3) be satisfied. Then the following 
statements are true: 

(1) The inverse problem (0.1)- (0. 4) has a solution in Qt; 

(2) The solution of inverse problem (0.1)-(0.4) is unique in Qt 01 where the 
number T (0 < T < T) is determined by the data of the problem. 
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Proof. Any solution of the equation (0.1) can be given by 

oo 

(2.1) u(x,t) = ^u k (t)X k (x), 

k=l 

where the functions Uk(t), k = 0, 1, 2, ... satisfy the following system of equations: 

u' (t) = F (t), 

u' 2k (t) + (2irkf a(t)u 2k (t) = F 2k (t), 



«2fc-i(*) + ( 2 ^ fc ) 2 o(*)«2fc-i(t) + ^ku 2k (t) = F 2 fc-i(t), fc = 1, 2, ... . 

Substuting the solution of this system of equations and initial condition (0.2) in 
(2.1), wc obtain the solution of the problem (0.1)-(0.3) in the following form 



u(x, t) 
(2.2) 



<Po + 



X (x) 



fe=i 

oo 

fe=l 

oo 



-(27rfc) 2 fa(s)ds f -(27rfc) 2 / a(s)ds 

^ 2fe e + / F 2k (r)e - dr 



^2fc(a0 



fe=i 



-(27T-fc) 2 fa(s)ds 

(f2k-i - 47T k<p 2k t) e 



Z 

J (F 2k _ 1 (T)-4nkF 2k (r)(t-r)), 



X 2 k-i(x) 

t 

(2-rrk) 2 f a(s)ds 



dT 



X 2 k-i(x). 



Under the conditions (1) of (A 2 ) and (1) of (A3) the series (2.2) and ~§x are 

fe=i 

uniformly convergent in Q T since their majorizing sums are absolutely convergent. 
Therefore their sums u(x,t) and u x (x,t) are continuous in Q T . In addition, the 

OO OO 2 

series Y and Yl ~§^z are uniformly convergent in Qt- Thus, we have u(x,t) G 
k=i k=i x 

C 2 ' 1 (Qt) H C 1 ' (Qt)- In addition, Ut(x,t) is continuous in Q T . Differentiating 
(0.4) under the condition (Ai) , we obtain 



(2.3) J u t (x, t)dx = E'(t), < t < T, 



and this yields 



(2.4) 



a(t) = P [a(t)} , 
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where 

2*b (*) + E ^F 2k (t)-E'(t) 



oo 

2 



(2.5) P[a(t)} = 



k=i 



-(2-Kk)* J a(s)ds t ( 27r fc) 2 / a(s)ds 

E 87rfc </> 2 fc e + / -F 2 fc(r)e - dr 

k=l \ 



Denote 



Co = 2 min F (t) + min (V -^-F 2k (t) ] - max £'(t), 
te[o,T] te[o,T] I 7rfc / te[o,T] 

Ci = 2 max F (t) + max [V -^-F 2k (t) | - min £'(*), 
te[o,T] te[o,T] V ^ 7rfc / *e[o,T] 

T oo T 

C 2 = J E'{t)dt + J^2fc - 2 J F (t)dt, 

k=1 

It is easy to verify that Cfc > 0, k = 1, 2, 3, 4 and C 2 < C 3 , such that 



T 

f ^2 f 2 / -(27rfc) 2 f o(s)ds f -(2nk) 2 ( a{s)ds \ 

j E'(t)dt+Y. j F °® dt ^ e ^ ^ e + y ^ 2fe(r)e H • 

' c=1 fc=1 V / 

Using the representation (2.4), the following estimate is true: 

o<^<«(t)<^. 

C3 G 2 

Introduce the set M as: 

M={ fl ( t )6C[0 1 3l:|<a(t)<|}. 
It is easy to see that 

P : M -> M. 

Show that the operator P is compact. Let Mi C M be an arbitrary bounded set. 
Since P(Mi) C M, then P(Mi) is uniformly bounded. Then, we have for a(i) G Mi 
and ii, i 2 G [0,T], 

( , 6) |fKllFfMl ,,||^ ^ MS - r .i 

where 

oo 2 

jf (*) = 2Pow + y; ^-p 2k {t) - ^(*), 
fc=i ™ 

/ -( 2 7rfc) 2 f a(s)ds /" -( 2 7rfc) 2 f a(s)ds \ 

N(t) = Y^8nkU 2k e i + j F 2k (r)e I dr\. 
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Using the estimate 



-(2-Kkf f a(s)ds -(2-nk) 2 f a(s)di 

e — e 



< (27rfc) 2 It! -t 2 \max.a(t), 

[0,T] 



we obtain 

(2.7) 

where 



\N(h) - N(t 2 )\ < 



(C 4 + C 5 ) ^- + C 6 



1*1 -*2|, 



oo „ oo / oo \ 

C 4 = ^4(27rfc) 3 ^ 2fe , C 6 = ^4(27rfc) 3 F 2fc (r)dr, C 6 = max ^87rfcP 2fc (i) 
fc=i £ fc=i te[0 ' T] \fe=i / 

To this end, take an arbitrary e > 0. 

Since if (t) is continuous in [0, T] , then 35 1 = 5i(e), V*i, i 2 € [0,T] (|*i -i 2 | < <*i) : 



(2.8) 
Let 



|A-(ti)-if(t 2 )|< 



C 2 £ 



(5 = min |#i(£), 



2((C 4 + C 5 )C 1 + C 2 C 6 )C 1 
From (2.7) for \t\ — t 2 \ < 5, we obtain 



(2.9) 



\N( tl )-N(t 2 )\<^e. 



Substituting (2.8) and (2.9) in (2.6) we get 

|P[o(ti)]-P[o(t 2 )]| <e- 

So, the set P{M\) is equicontinuous. Then P(M\) is a compact set and the 
operator P is compact and maps the set M onto itself. Employing the Schauder's 
Fixed Point Theorem, we have a solution a(t) e C [0, T] of the equation (2.4). 

Now let us show that there exists Qt (0 < To < T) and the solution (a, it) 
of the problem (0.1)-(0.4) is unique in Qt - Suppose that (b, v) is also a solution 
pair of the problem (0.1)-(0.4). Then from the representation (2.2) and (2.4) of the 
solution, we have 



, -(27rfc) 2 fa(s)ds -(27rfc) 2 fb(s)ds 

(2.10) u(x,t)-v(x,t) =2^ 2fe ( e o _ e o \ X 2k (x) 

fe=i 



-(27rfe) 2 fa(s)ds -(Ink) 2 f b(s)ds 

( r) | i - e I '/r A-j/, (./•) 



, -(27rfe) 2 fa(s)ds -( 2 Trkf f b(s)ds 

2_^VP2k-i - ^kip 2k t) | e o -e ° | X 2fc -i(aO 

fc=i 
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, , , -(27rfe) 2 fa(s)ds -(2irk) 2 fb(s)ds 

+ E / (^fe-i(r) - 4nkF 2k (r)(t - r)) e ' - e ' | dr | X 2fe 



(2.11) 
where 



a(t)-b(t) = P[a(t)]-P[b{t)} 



P[a(t)}-P[b(t)] 



k=l 



-(2vkr f a(s)ds t -(2irk) 2 f a(s)ds 

£ 87rfc U> 2fe e o +/F 2fc (T)e - 

fe=i \ o 



2f b(t) + E £*2 fc (t) - 



fc=i 



, -(2 OT fc)2/6( s ) ds t -(2-KkY fb(s)ds 

E 8rrfc ^ 2fc e ° + /F 2 k(T)e - 

fe=i \ o 



The following estimate is true: 

/ oo 

(<2F (t)+j:±F 2k (t) + \E'(t)\ 



\P[a(t)]-P[b(t)]\ < 



k=i 



CI 



( 



E 87rfcv5 2fe 



fc=i 

oo T 



-(27rfc) 2 fa(s)ds -(2irk) 2 fb(s)ds 

e — e 



E 87rfc/^ 2fe (r) 

fc=l o 



-(27rfe) 2 /a(s)ds -(27rfc) 2 fb(s)ds 

e T — e T 



aV 



Using the estimates 



-(2nk) 2 J a(s)ds -(2nk) 2 j b(s)ds 

e — e ° 



-(2-Kkf f a(s)ds -(2ixkf f b(s)ds 

e t — e 



< (27rfc) 2 T max \a(t) - b(t)\ , 

0<t<T 



< (27rfc) 2 T max \a(t) - b(t)\ , 

0<t<T 



we obtain 



max IP [ami - P \b(t)}\ < a max a(t) - . 

0<t<T 0<t<T 

Let a 6 (0, 1) be arbitrary fixed number. Fix a number To, < To < T, such that 

Ci(<7 4 + C 6 ), 



<3 



-To < «• 
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Then from the equality (2.11) we obtain 

ll a " 6 llc[0,T ] < a \\ a - b \\c[0,T ]> 

which implies that a = b. By substituting a = b in (2.10), we have u = v. □ 
Theorem has been proved. 

3. Continuous Dependence of (a, u) upon the data 

Theorem 2. Under assumption (Ai) — (A3), the solution (a,u) of the problem 
(O.l)-(O.S) depends continuously upon the data for small T . 

Proof. Let $ = {(p, F, E} and $ = {jp, F, E} be two sets of the data, which 
satisfy the conditions (Ai) — (A 3 ) . Let us denote ||$|| = C 1 1 V 1 1 c?-* [0 1] + ll^llc 4 .°(Q T ) + 
||-E|| C i[ T j). Suppose that there exist positive constants Mj, i = 1, 2, 3 such that 

IMIc*[o,i] < M L 11^11^.0(0^) < M 2, PIIcmct] < M 3 , 

\\n c *io,i] ^ M i. FL«.o(<jd < M 2, FIIcmct] ^ M 3- 

Let (a, u) and (a,u) be solutions of the inverse problem (0.1)-(0.4) corresponding 
to the data 3> and <I>, respectively. According to (2.4), 



a(t) 



2*b (*) + E ^E^fcW - S'(t) 
fc=i 



-(2^) 2 /a( s )d s * -(27rfe) 2 /a( 5 )d S 

E Stt/c ip 2k e ° + j F 2k { T )e - dr 

fe=i \ 



a(t) 



2^o(i)+ E 4^ fe (t)-SW 
fc=l 



-(27rfe) 2 fa(s)ds t_ -(2-xk) 2 fa(s)ds 

E 87r/c ^ 2fc e +/F 2fe (r)e - dr 

fc=i \ 



First let us estimate the difference a — a. It is easy to compute that 



~ -(27Tk) 2 fa(s)ds _ ~ -(2TTk) 2 fa(s)ds 

F (t)}^8Trklp 2k e ° - F (t)}^8TTk(f 2k e 



k=l 



fe=l 



< M 4 ||a-a|| +M 5 \\<p-W\\ + M 6 \\F - F\\ 



~ f_ -{2-nkf }a(s)ds _ ~ f -(27rfc) 2 fa(s)ds 

F (t)J2^k F 2k (r)e ' dr - F (t) ^ SttA; / F 2k (r)e ' dr 



k=i 



< M 7 \\a-a\\ +2TM 5 \\F - F\\ 

' " "C[0,T] °ll H C4.0(Q T ) 



fe=l 
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2 _ -(27rfe) 2 fa(s)ds ^2, 2 — , -(27rfc) 2 fa(s)ds 



^2—F 2k {t)^^kip 2k e 



fe=i 



fe=l 



fe=l 



fe=l 



^4 „ ^5 „ ,, ^6 n —II 



oo oo t -(27rfc) 2 f a(s)ds 

E ^F2k(t) E 8nk J F 2k (r)e t dr 

fe=i fe=i o 



- (27rfc) 2 fa(s)ds 



k=l 



< 



T 2 M 7 



k=l 

TM 5 



..a — a + 

6 " m c[o,t] 3 



F - F\ 



C 4 .0(Q T ) ' 



~ -(2nk) 2 fa(s)ds _, ~ -(27rfe) 2 /a( S )d S 

£'(i)^87r^ 2fe e - - £ (t) }J 87rfc^ 2fc e 

fc=i fe=i 

< M 8 ||o-a|| + M g \\ip-Tp\\ + M 6 ||£ - Ell 

- 011 "c[0,T] J ^ M 04[0,1] ° 11 N Ol[0,T] 

t t t 

^2, /•_ -(27rfc) 2 / a(s)ds _, ~ /" 

£'(t) J] STrfc J F 2k {r)e ' dr-E(t)J2 J F 2k (r) 



k=l 



(2tt/c) 2 fa(s)ds 

e dr 



k=l 



< M w \\a-a\\ + TMq \\F — F\\ + TM 5 \\E — E\\ 

where M k , k = 4, 12, are some constants. 

If we consider these estimates in a — a, we obtain 



(1-Mn) \\a - all < M12 LE - E\\ +\\w-W\\ + \\F — F\\ 

^ n ' M C[0,T] - \JI Hcl[0,T] ^"c4[0,l] 11 H C4,0(Q T ) 

The inequality M\\ < 1 holds for small T. Finally, we obtain 



— ,, Ml 9 

la — all < Mi 3 , Mi 3 = ^— . 

" "c[0,T] - II II ' (1 - Mil) 

From (2.2), the similar estimate is also obtained for the difference ti-uas 



\c{Q T ) 



'14 



□ 



4. Numerical method and an example 

We consider an example of numerical solution of the inverse problem (0.1)- 
(0.4). We use the finite difference method with a predictor-corrector type approach, 
that is suggested in [2]. Apply this method to the problem (0.1)-(0.4). 

We subdivide the intervals [0, 1] and [0,T] into M and N subintervals of equal 
lengths h = and r = ^ respectively. Then we add a line x = (M + 1) h to 
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generate the fictitious point needed for dealing with the second boundary condition. 
We choose the Crank-Nicolson scheme, which is absolutely stable and has a second 
order accuracy in both h and r. ([9]) The Crank-Nicolson scheme for (0.1)-(0.4) is 
as follows: 



1 I j'+i 

— \U; — U 



(4.1) 



F: 



(4.2) 



«° = <t>i, 



(4.3) 



u n = u 



Mi 



'M-l 



= U 



M+l' 



(4.4) 

where 1 < i < M and < j < N are the indices for the spatial and time steps 
respectively, u\ — u(xi,tj), <^> i — f{xi), F- = F(xi,tj), xi = ih, tj = jr. At the 
t = level, adjustment should be made according to the initial condition and the 
compatability requirements. 

(4.1)-(4.4) problem forms M x M linear system of equations 



(4.5) 
where 



AU 1+1 = b, 



u 3 - ">;) , i < 3 <n, b = (b u b 2 ,...,b M y j 



A = 



-2(1 + R) 1 

1 -2(1 + R) 1 ... 

1 -2(1 + R) 1 ... 



R 



2h* 







1 -2(1 + R) 
2 



r(oi+ 1 +ai)' 3 ~ 0, 1, N, 

h = 2(1 - R)u{ - u{ - u j M - Rt(F( +1 + F(), j - 0, 1, N, 
bu 
h = 



-2?; J 

AU M-1 



- 2(1 - R)v? M - RriF^ 1 + F M ), j = 0, 1, N, 



-2{l + R) 



_ 1 +2(1-R)i 



! - Rr(Fl +l +F?),i = 2, 3, M-l, j = 0, 1, N. 



Now, let us construct the predicting-correcting mechanism. First, integrating 
the equation (0.1) respect to x from to 1 and using (0.3), (0.4), we obtain 



(4.6) 



a(t) 



-E'(t) + F(x,t)dx 



u*(0,t) 

The finite difference approximation of (4.6) is 
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(Et) 3 + (Fin)*') h 



u 



where (Et) J = E'(tj), {Fin) j = J* F(x, t 3 )dx, j = 0, 1, N. For j = 0, 

_ (-(Et)° + (Fin)°)h 

and the values of <j> i provide us to start our computation. We denote the values of 
o J , u{ at the s-th iteration step a^ s \ v^ s \ respectively. In numerical computation, 

since the time step is very small, we can take a^ +1 ^ = a 3 , v^ +1 ^ = uj, j — 
0,1,2, ....N, i — 1,2,..., M. At each (s + l)-th iteration step we first determine 
fl j+i(s+i) f rom the formula 



+ (- (Ety +1 + {Finy +1 ^ h 
Then from (4.1)-(4.3) we obtain 



,/+!(«)_ „i+i(«) 

"1 "0 



(4.7) + (t^^ - 2«f^ +1 > + + 1 + Ff) 

(4-8) ^ +1(s) =< +1(s) , 
(4.9) «^- ( i a) =«^ ( i S) . « = 0,1,2,.... 

The problem (4.7)-(4.9) can be solved by the Gauss elimination method and u^ +1 ^ 8+1 ^ 
is determined. If the difference of values on two iteration reaches the prescribed tol- 
erence, the iteration is stopped and we accept the corresponding values a i+1 ^ s+1 \ 
ul +1(s+1 \i = 1,2,..., M) as a j+1 , ul +1 (i = 1,2,..., M), on the (j + l)-th time step, 
respectively. In virtue of this iteration, we can move from level j to level j + 1. 

Example. Consider the inverse problem (0.1)-(0.4), with 

F(x,t) = ^— cxp(— t) + 47rcxp(3t)^ cos27rx + (27r) 2 (l — x) sin27rxexp(3i), 

1 1 

ip{x) = (1 - x)sin27nc, E{t) = — exp(-i), T = -. 

2ir 4 

It is easy to check that the exact solution is 

{a(t), u(x, t)} = < 2 + ex P (4^) ! (1 — x) sin 27t.t cxp(— t) > . 

1 ( 27r ) J 

We use the Crank-Nicolson scheme and the iteration which are explained above. 
In result, we obtain Table 1 and Table 2 for exact and approximate values of a{t) 
and u(x, t). The step sizes are h — 0.005 and r = \. 
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Table 1. 


The some values of a(t) 




Exact 


Approximate 


Error 


Relative Error 


1.0354 


1.0261 


0.0093 


0.009 


1.3685 


1.3276 


0.0409 


0.0299 


1.6658 


1.6308 


0.035 


0.021 


1.9698 


1.9487 


0.0211 


0.0107 


2.0798 


2.0644 


0.0154 


0.0074 


2.2068 


2.2008 


0.006 


0.0027 


2.3186 


2.322 


0.0034 


0.0014 


2.4727 


2.4915 


0.0188 


0.0076 


2.5981 


2.6304 


0.0323 


0.0124 


2.7301 


2.7778 


0.0477 


0.0175 
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Table 2. The some values of u(x, t) for T = 70 



Exact 


Approximate 


Error 


Relative Error 


0.3633 


0.3486 


0.0147 


0.0347 


0.5952 


0.5933 


0.0019 


0.0044 


0.6506 


0.6523 


0.0017 


0.0001 


0.4054 


0.4084 


0.003 


0.0017 


0.3402 


0.3415 


0.0013 


0.0016 


0.2563 


0.2554 


0.0009 


0.0085 


-0.2157 


-0.2298 


0.0141 


0.0457 


-0.1430 


-0.1517 


0.0087 


0.0066 


-0.1176 


-0.1252 


0.0076 


0.0019 


-0.0754 


-0.0812 


0.0059 


0.0276 
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